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ABSTRACT 



In this paper we specifically studied one aspect of foreground primordial mat- 
ter density perturbations: the relative gravitational time delay between a pair of 
light paths converging towards an observer and originating from two points on the 
last scattering surface separated by the physical scale of an acoustic oscillation. 
It is found that time delay biases the size of acoustic oscillations systematically 
towards smaller angles, or larger harmonic numbers £, i.e. the mean geometry 
as revealed by CMB light becomes that of an open Universe if f2 = 1. Since the 
effect is second order, its standard deviation 51/ i ~ (5$)^ where (5$)^ ~ 10^^ 
is the normalization of the primordial matter spectrum P{k), the consequence is 
too numerically feeble to warrant a re-interpretion of WMAP data. If, however, 
this normalization were increased to 5$ > 0.01 which is still well within the per- 
turbation limit, the shift in the positions of the acoustic peaks would have been 
substantial enough to implicate inflationary ACDM cosmology. Thus fl is not 
the only parameter (and by deduction inflation cannot be the only mechanism) 
of relevance to the understanding of observed large scale geometry. The physics 
that explains why 5$ is so small also plays a crucial role, but since this is a 
separate issue independent of inflation, might it be less artificial to look for an 
alternative solution to the flatness problem altogether? 

Subject headings: 



1. Introduction 

Among the so-called 'secondary' physical mechanisms that re-process the CMB anisotropy, 
such as gravitational lensing, time delay, and the Sunyaev-Zel'dovich effect, time delay by 
foreground inhomogeneities in the matter distribution appears to be investigated least. In 
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the recent period a detailed treatment of the problem was provided by Hu & Cooray (2001, 
hereafter HCOl), although the general framework for calculating gravitational perturbation 
effects on the CMB as published in the review of CMB lensing by Lewis & Challinor (2006) 
could also be employed to carry the study further. 

In HCOl the authors found an infra-red logarithmic divergence in the variance of the 

absolute (or total) time delay along a randomly chosen direction to the LSS. They 'renormal- 
ized' this infinity by subtracting the contribution from the 'monopole' term of the matter 
power spectrum, corresponding to the removal of a constant uniformly across the sky. Nev- 
ertheless, the remaining (finite) quantity still carries the divergence trend. More precisely 
the variance is dominated by long wavelength fluctuations, prompting HCOl to consider it 
as an effect of very large coherence length, A£ fa 2 which does not affect our interpretation 
of the CMB anisotropy, because the net outcome is simply a gentle and arbitrary distortion 
of the spherical shape of the LSS, completely neghgible over the size of one (or a few) cycles 
of CMB acoustic oscillations. 

The coherence length inferred by HCOl should be viewed with some caution, however, 
because the large delay excursion of ~ 1 Mpc/c calculated there, which involved only the 
zeroth order term, the path integral of the perturbing potential itself, stems from the part of 

the matter power spectrum which carries the scale-invariant Harrison-Zel'dovich dependence 
P{k) ~ k, i.e. there is a danger that the large coherence length may simply be due to the 
infra-red divergence of the variance rather than any genuine physical scale in the matter 
spectrum. To find the coherence length of relevance to the question of CMB anisotropy 
distortion, it is necessary to pursue the perturbation expansion to the next two orders. Not 
only are both results free from divergences, but also it is only through a comparison of these 
two terms that the true coherence scale for variations in the relative delay between two 
hght paths separated by a small angle 9 would become transparent. We shall find that this 
scale is defined by a physically significant parameter, viz. the characteristic wavenumber at 
which departures of P{k) from the Harrison-Zel'dovich behavior occurs for the first time. 
The consequence is that appreciable distortion of the LSS radius, with both amplitudes and 
wavelengths on par with the dimension of the primary acoustic oscillations at the time of 
last scattering, can exist in principle. 

2. Perturbation in the gravitational potential from the 2dFGRS/WMAPl 

power spectrum of primordial matter 

Although in the k ^ limit the matter power spectrum has the form P{k) ~ k, the 
behavior of P{k) at large k is more complicated than an exponential cutoff. It is possible. 
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however, to break down any general P{k) into constituent terms, each of the form aike"^"^^, 
and sum up the time delay fluctuation contributions from all the terms, because P{k) has 
the meaning of a variance, i.e. it too is additive. We may therefore write 

P{k) = A{a^ke-^'^ + a2ke-^^^ + •••), with J] = 1. (1) 

i 

This empirical representation of P{k) is not the same as the more commonly used ones 
(e.g. Efstathiou, Bond, and White 1992) but, as shall be seen in section 3, the exponential 
form reveals coherent length scales of foreground effects in a transparent way; in any case, 
provided our formula for P{k) fits the observational data (see below) the detailed structure 
of the terms used to model the spectrum is of no significance. 

The resulting value of A in Eq. (1) that we shall obtain is 

A = 3.276 X 10^ Mpc^ (2) 

and leads, by Eq. (A-6), to 

5$ 3 X 10-^ (3) 

for a = 0.3, I^a = 0.7, and h = 0.7 cosmology (Bennett et al 2003, Spergel et al 2007). 
This agrees well with the CMB temperature modulation of 3ST/T at small k as measured 
by WMAP (see e.g. Bennett et al 2003), as it ought to, because from Eq. (A-5) 

i^^) =^^.id: = }'^l{^] ' (4) 



k^Q dink k^o \ Ti 

where the final step is explained in the material around Eq. (18.14) of Peacock (1999). The 
consistency between 5$ as derived from our z — Q matter spectrum and the large scale 
CMB anisotropy re-assures us that any corrections we ignored, such as the effect of vacuum 
domination at z < 0.3, are indeed minor. 

If P{k) has the simple form involving only the first term of Eq. (1) with oi = 1, we may 
work out from Eq. (A-5) the correlation function for the pcrtiubing Newtonian potential $ 
that arises from the linear growth of primordial density contrasts, as 

($(r)$(r')) = j '^e^k<^-^')p{k) = j ^^Mr-r')^-bk ^ 

where $ is assumed to be time-independent and Gaussian distributed, and in the final step 
use was made of Eqs. (A-4) and (A-6). If P{k) is given by the full Eq. (1) instead, then the 
rightmost side of Eq. (5) will be a sum of similar terms, each carrying the exponent hi and 
with {5^Y replaced by aj(5$)^. 
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Implementing now the observed power spectrum, the most up to date data are from 
the 2dFGRS galaxy survey (Cole et al 2005) after they are deconvolved and aligned with 
the WMAPl normalization by setting the (t| parameter to (jg = 0.74 (Sanchez et al 2006). 
Wc found that the resulting dataset can adequately be fitted with a function for P{k) of 
the form given by Eq. (1) and involving three exponential terms, with the value of A as 
already quoted in Eq. (2) and the values of and bi {i = 1,2,3) as shown in Table 1. This 
best-fit spectrum, which closely follows that of WMAPl's ACDM model (Spergel et al 2003) 
is plotted in Fig. 1. 
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Fig. 1. — The best P{k) model of the WMAPl normahzed 2dFGRS data as given by Eq. 
(1) with A = 1 and the remaining parameters as shown in Table 1. The Hubble constant 
assumed is h = 0.7. 
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3. Cosmological time delay veiriation from primordial density contrasts 

We proceed towards calculating the excursion in time delay along two light paths of 
equal lengths. Our starting point is a flat Universe having its average density at the critical 
value, as required by the WMAP observations (Bennett et al 2003, Spergel et al 2007). We 
employ Cartesian comoving coordinates and the conformal time coordinate rj, so that in a 
perturbed flat FRW space light propagates along the null geodesies of the metric 

ds^ = (1 + 2^)d7]'^ - (1 - 2^){dx'^ + dy^ + dz^), (6) 

as though the expansion factor a{r]) plays no role. Let light signals arrive at the observer's 
origin from a source at (x, 0, 0). Although the signal is coming towards us, we can, by optical 
reciprocity, solve the equation starting at x = and following its path backwards. 
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Table 1: Parameters (with l-cr errors) for the best P{k) model of the WMAPl normahzed 

2dFGRS data {h = 0.7). The form of the model is given by Eq. (1). Three exponentials 

were needed to fit the data. 
Oj Gi error bi (Mpc) bi error (Mpc) 

0.883 017 128.0 11.0 

0.114 0.02 39.3 3.80 

0.0017 0.001 9.01 0.80 
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Suppose a signal arrives from a direction making a small angle w.r.t. the x axis. This 
may correspond to some 'off- axis' point on the LSS in the case of the CMB. The conformal 
time of travel for our light signal is, from Eq. (6), 



V 



with an ensuing time delay of 



[1 - 2$(a;', dx')]dx\ 



f 2$(x', ex')dx', 
Jo 



(7) 



(8) 



where x and x' can be present day physical distances if we set a{z = 0) = ao = 1 (here we 
ignored the geometric time delay, with is comparatively negligible as pointed out by HCOl). 
Moreover, if we use the current value of $, as in the previous section, to calculate r, then 
6t = 5t/a{t) = St. The 'absolute' variance of the time delay along a random direction to 
the LSS 

{t\0))= r 2dx' r 2dx"mx',0x')^x",0x")), (9) 
Jo Jo 

is logarithmically divergent because the quantity ($(r)$(r')) is, by Eq. (5), of the form 
dk'/k' at small k. As stated in the beginning of the paper, HCOl renormalized (t^(0)) 
by subtracting an infinite constant from it. 

Of more relevance to understanding the CMB acoustic peaks is the relative delay in the 
arrival time between the above signal and another light signal emitted simultaneously from 
the same distance, but along the 'on-axis' direction — 0, i.e. 



t{0)-t{O)= 2x'0 ■V^{x',0)dx' + x'\0 ■Vy^{x',0)dx' + 
Jo Jo 



(10) 



where V is the gradient operator transverse to the vector x. assuming for the time being 
that the form of P{k) is given by Eq. (A-4), we can construct the correlation function with 
the help of Eq. (5), as 



(5$)2 f d^k,, jk-r„-bk 



^^-kikjc'" ■ e 



b /r 
1 arctan - 



+ 



3b 

— arctan 
r \b 



\bJ r^ + b'^j) 



(11) 



where 



r' - r", 



(12) 
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and the indices i,j denote two orthogonal components in directions transverse to x. This 
enables us to derive the lowest order term for the variance in the relative time delay [(^t(^)]^ = 
{[r{6) — t(0)]^), viz. the first term on the right side of Eq. (10), as 



{S^ fd^ I 2x'dx' I 2x"dx" 
'o Jo 



lb /r 
— arctan - 



where r = \x' — x"\, consistent with Eq. (12). 

Transforming now to the new variables x = {x' + x")/2 and x — x' — x", the resulting 
integrand is symmetric in x, so we can restrict the range of x to positive values, introducing 
an extra factor of 2. Thus 



M9)f 



2{6^yd^ / dx 
Jo 
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(13) 



By restricting ourselves to the limit x ^ b ((appropriate to emission distances x ^ 6i ~ 0.1 
Mpc - see Table 1 - the CMB LSS clearly satisfies this criterion) , it becomes straightforward 
to complete the calculation, because only one term stands out. The result is 



MO)? - ^(<^$)^0^y arctan (^) = fm'e'j. 



(14) 



Note that according to Eq. (14) the variation in the time delay difference between two points 
A and B on the LSS subtending an angle 6 at the observer O is St{9) ~ 9. This behavior 
indicates that we are in the regime of coherent delay, i.e. provided 9 is sufficiently small the 
two rays sampled a primordial matter potential gradient which may be regarded as constant. 
Thus, if in Figure 2a a third point S3 on the LSS is coUinear with two other points Si and S2, 
then a constant gradient would imply equality between the 'Si to S2' and 'S2 to S3' relative 
delays, i.e. both are St{9). Consequently the net Si to S3 delay will be 26t{9). Since the 
angular separation between Si and S3 is 26 at O, we have 5t{29) = 25t{9), fully consistent 
with the 5t{9) ~ 9 dependence. 



Fig. 2. — Illustrating the phenomenon of coherent time delay. 2a: simultaneously emitted 
light signals from three symmetrically located and equally distant sources would have arrived 
at points A, B, and C also at the same instance if there were no gravitational time delay. 
Provided d is small enough that only the first spatial derivative of the potential is responsible 
for any delay, then the signals from Si and S2 will e.g. have reached A' and B' when that 
from S3 is at C'=C, where AA' = 2BB'. The result is a net displacement of the angular 
position of all three sources by the same amount to the right. 2b: If we reverse the light 
paths by letting the observer be a single source and the former (Si, S2, S3) be replaced by 
(B, O, A), then the source will appear shifted in the opposite direction (i.e. to the left) w.r.t. 
the observer whose telescope aperture is AB, as it should do. 
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Turning to the term, which takes into account the difference in the potential 

gradient between the two rays, its exact mathematical form can be derived directly from 
the last integral of Eq. (10), but in this paper we will present instead another method of 
calculation which also provides unusual insights to the relationship between time delay and 
gravitational lensing. Such material is to be found in sections 4 and 5. For now we simply 
give the result, i.e. an expansion of [(5r(^)]^ to include the term. It is 

Me)r - Y^m'o'^ (i - for e < e^- (is) 

or, if P{k) has the more general form of Eq. (1), the contribution from the ith term to the 
variance (note the label i here has a different meaning from that in Eq. (11)) will be 

[6r.iO)r - (l - for e < ei, (16) 

where in Eqs. (15) and (16) the parameters 9m and 6^^ are yet to be defined below - 
see Eqs. (17) and (20). From Eq. (15) we see that the higher order [6^) term, once it 
assumes importance, will halt the linear rise of St{6) with 6. Eventually, when 9 becomes 
large enough, all the higher order terms of Eq. (10) will take their place to ensure that 
[5r(6')]^ reaches constanc}|l], as it must do, because two widely separated rays are uncorrelated. 
Thus, once the 6'^ term is no longer negligible, incoherence takes over. A reasonable (and 
conservative) way of estimating the maximum (or 'plateau') value of [(5r(6')]^ is to find the 
angle at which the first ^-derivative of the right side of Eq. (15) vanishes. This occurs when 

= VlO-, (17) 

X 

at which point 6t reaches the constant (saturation) value of 

[5r{0)]' = Mej]' = ^(5$)'x6 for 6 > (18) 

Here 9m may be defined as the coherence angle for time delay. It is also closely related to 
the expansion parameter of the perturbation expansion of Eq. (15), because the ratio of the 
second order to the first order term of Eq. (15) is 26m- When a more general form of P{k), 
like Eq. (1), is considered, we have 

[Srm? = [5r.iei)r = ^m'^aA for 9 > 9^, (19) 



^Strictly speaking there may be a gentle rise with 9 because of the logarithmic divergence problem 
mentioned in section 1. This has no significant impact on any of the results presented in this paper, however. 
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where 

C = ^ (20) 

for each component i of P{k) in Eq. (1). The total variance of time delay is then the sum 
over i of all the individual variances [5Ti{d)Y as given by Eqs. (16) and (19). Obviously, for 
rays separated by angles 9 > 9m the time delay is always completely incoherent (or random) . 
As will be demonstrated in the next section, the significance of the transition angle 9m is 
that for separation 9 < 9m, the variance for the coherent (or systematic) and incoherent 
delay is given by the first and second order terms of Eq. (15) respectively. 
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Fig. 3. — 1-0" relative time delay 5t^{9) between two rays from the LSS due to the various 
terms i of primordial density perturbation of Eq. (1). In the limit 9 < 9]^^ where 9]^ given 
by Eq. (20), 5t^{9) is given by Eq. (16). When 9 > 9\^, 6t\9) saturates to the value in Eq. 
(19). Long dashes correspond to the first term i = 1, and the dotted line to the second term 
i = 2 (the i = 3 term was ignored due to its insignificance). Values of and bi are shown in 
Table 1. 
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How does time delay affect CMB observations? Wfien ligfit signals emitted during 
the same redshift z = zlss from two points on the LSS which subtend the angle 9 at the 
observer are to arrive simultaneously and yet the times of flight are delayed w.r.t. each 
other, the signals would necessarily have covered different distances - the slower one must 
have undertaken a shorter journey. As illustrated in Figure 4, the consequence is a 'tilt' of 
the LSS in that vicinity. Hence, an angle 6 which normally corresponds to an anisotropy 
on the LSS spanning the comoving distance xd would now be randomly mapped to a larger 
distance x{9 + 59) where, from Figure 4, 



with the contribution to {5TiY from the ith term of Eq. (1) given by Eqs. (16) and (19). 
The net outcome is a systematic shift of the CMB acoustic peaks, which are features with 
fixed comoving scales, towards smaller angles or higher harmonic numbers I. 

It should be mentioned that another possible test of cosmological time delay concerns 
two light paths propagating through different parts of the Universe, but connecting the same 
source with the observer, i.e. a strong gravitational lensing scenario, under which multiple 

images occur. If this source undergoes flaring behavior, a time lag between the light curves 
of a pair of multiple images might result from perturbations in the distribution of primordial 
matter. In practice, however, this is a sensitive test only if the lowest order {5t{9) ~ 9) term 
contributes to the delay. As it turns out, this term does noi play a role because the potential 
time delay we have hitherto been considering is cancelled by another effect - the geometric 
delay - which is important under this scenario (see Bar-Kana 1996 and Seljak 1994). Since 
the higher order terms are too feeble to be measured, no meaningful constraints can be 
provided here. 




(21) 
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Fig. 4. — The effect of CMB time delay on the appearance of the acoustic peak hot and cold 
spots. Two proximity locations on the observed LSS are the points A and B. Were time 
delay absent, they would have been the points A and C. Here the distance EC is typically 
5t{9), the quadrature sum (over all i) of the STi{6) contributions as given by Eqs. (16) and 
(19), because it is the relative delay between the light signals emitted at A and C, with 
the path AO suffering from more delay than BO in this case. Since any LSS plane samples 
the same distribution of hot and cold spots as any other, the tilting of the LSS by time 
delay invariably means an observer charted slightly larger physical distance on the LSS for a 
certain angular separation between two temperature sensors. This systematically causes all 
acoustic peaks to move towards smaller 6, or larger i, in the manner of Eq. (21). 
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4. Deflection angle fluctuation from primordial density contraists 

A significant gain in our general understanding of foreground effects is afforded by 
examining also the question of CMB lensing. In addition to the x-axis of section 3, which 
points along a radial direction from the observer's origin back to the LSS, we now introduce 
two more Cartesian (comoving) coordinates y, both measuring distances transverse to the 
light path. Tracing the path towards the LSS as before, the equation of motion for a hght 
signal is 

= -2V*. (22) 

Let a signal arrive from some direction which makes a (small) angle d with the x axis. Then 
the solution to (28) is 

y^dx + dy^ex-2 f dx' (x - x')V^x' ,ex'), (23) 

Jo 

with appropriate initial conditions. 

From Eq. (30) emerges that the correlation between the values of the deviation angle 
6y/x for different values of 0, 



x"^ 
4 

x"^ 



r dx'{x - x') r dx"{x - x"){V[^x\ d'x')V;^x", e"x")). (24) 
Jo 

By means of Eqs. (24) and (11), one can expand the ensuing correlation function C{d) as a 
Taylor series (see Appendix B for details on intermediate steps), i.e. 

^ cm . Mi^MzM . (i - . . . .) . (25) 

where the sum represented by a repeated index is over the two transverse dimensions. 



5. Reciprocity of light propagation: the relationship between time delay and 

deflection 

An interesting connection exists between sections 3 and 4, in that the key formulae 
derived in each section, viz. Eqs. (15) and (25), are closely related to each other because of 
a physical reason - the reciprocity of light propagation. 

We first consider the scenario of a source S emitting light signals that enter an observer's 
circular telescope aperture at the extremeties A and B such that the angle ASB is 29, see 
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Fig 2b. Let us assume that when the Universe is homogeneous, the wavefronts are parallel 
to the aperture as they arrive, i.e. the signals reach points A and B simultaneously. We 
further suppose that the same is true for two points C and D on the extremeties of another 
aperture diameter perpendicular to the AB line - the source is 'on-axis'. 

In the presence of primordial matter density perturbation, let the signal propagating in 
the SB direction reach the point B' when the SA signal has already arrived at A, i.e. the 
SA signal suffered from a smaller delay. The outcome is a tilting of the wavefronts as they 
enter the aperture - the source is now seen to have moved 'off-axis' along the AB line to 
a new position leftwards of the center O. The angle of tilt is 59 = 2Sr{9)/d in the limit of 
coherent time delay (i.e. small 6) where St{6) is the distance BB' and d is the diameter AB. 
Repeating our above argument to the points C and D, and noting that the relative delay 
between the SC and SD directions is independent of that between SA and SB, we realize 
that the variance 

iSe? ^ M (26) 

applies to the overall shift in the position of the source on the two dimensional sky. Of 
particular interest is the fact that for a given source distance x the quantity 

e = — ^d. (27) 

Thus, in order for the positional shift to be the same amount irrespective of aperture size d 
we must have 5r{6) ~ 9 where 9^ — 6-0. From section 3 we saw that this condition holds 
only when 9, hence ti, is small. 

To link sections 3 and 4, however, we have to consider a second scenario, under which the 
distribution of primordial matter remains the same as before, but A, O, and B are now three 
simultaneously emitting sources and S is the observer (equipped with a small telescope). 
This reversal of the light paths converts Figure 2b into Figure 2a. By the reciprocity of light 
propagation, the signals must arrive at S in such a way that the change in the positions of 
the three sources as perceived by S after the 'turning on' of the matter perturbation involves 
an angular shift by the same amount —56 for each source. Yet according to section 4, this 
shift has a variance given precisely by the correlation function C{0), where 6 is the angle 
the two sources subtend at the observer S, viz. 

.^Mi«M|!(zM = cW. (28) 
with C{0) being given by Eq. (25). 

If we now take the limit ^ ^ 0. By Eqs (28) and (22), the latter in section 4, we have 

(6ef = <^ = '-^m%, (29) 
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which is a constant shift (of O w.r.t. A, and B w.r.t. O) independent of d and 9. Such a 
behavior is also consistent with the requirement stated at the end of the previous paragraph 
when we considered our first scenario. Thus, from Eqs. (26), (27), and (29) we deduce that 

Me)r = fm'e'j, (30) 

which is the low 9 limit of the variance in the relative time delay between two light paths 
separated by angle 0, as derived in Eq. (15) of section 3. Thus, it is now clear that sections 
3 and 4 can be unified by the principle of light reciprocity. 

Of even more interest, however, is the regime of larger 9, where incoherence between 
the two rays becomes important. Here we already saw from Eq. (25) of section 4 that C(0), 
hence {SdY by Eq. (28), is no longer a constant, but decreases away from constancy as 
9 increases. This decrease is expected, because when 9 becomes sufficiently large the two 
light paths are independently perturbed by different primordial density fluctuations, i.e. the 
correlation function {6yi{^0)6yi{—^0)) , hence {60)^ by Eq. (28), assumes relatively small 
values. The two sources shift their positions randomly w.r.t. each other, leaving behind 
little net systematic drift of their centroid across the sky (according to S). In fact, when 
we compare Eq. (25) in full with Eqs. (27) and (28), we obtain an expression for [St{9)]'^ 
in complete agreement with Eq. (15). This consistency provided an important cross-check 
of the robustness of our analysis, and explains why the C2 term in the variance represents 
incoherence effects - it is precisely this term that decorrelates C{0). To labor upon this point 
even further, we observe that the variance of the random relative deflection between the two 
rays: 

- Sy{-m') = 2[C(0) - C{0)], 

has the C2 coefficient in is leading term. Thus, it is completely clear that while the Cq 
coefficient concerns absolute deflection and coherent delay, the physics of C2 is relative 
deflection and incoherent delay. 

A word of caution, however, before we leave this section. Although the phenomena of 
time delay and deflection are two sides of the same coin, the 'flipping of the coin' involves 
reversing light path arrows, i.e. for a given set of arrows we should not conclude that delay 
presents no new physical effects than those ensuing from lensing. In particular, for the CMB 
which is emitted from a three dimensional distribution of sources, a varying time delay on 
simultaneously observed signals from different directions allows depth to play a role in the 
problem - as is already explained towards the end of section 3. This depth effect cannot be 
reproduced in any way by lensing. 
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6. Time delay distortion of the CMB acoustic peciks 

We finally return to the original subject of this paper, viz. the degree to which time 
delay by foreground primordial matter re-processes the CMB primary anisotropy. In Eq. (1) 
and Table 1, we consider only the first two terms in our breakdown of density fluctuations, 
viz. i = 1, 2 of Eq. (1). Their coherence lengths are, by Eq. (20), e]^ = ei^ 0.0289 or = 
108.8, and = ^2 = 0.00886 or £2 = 7r/^2 = 354.4. Using Eq. (21), then, one derives the 
gaussian smoothing width of CMB anisotropy power, due to the i — 1 perturbation only, as 

- 4.5 X 10- X (5$/3 X 10-T X I (31) 



e ' ' ' ' \ 2(1 - 600^2), e<ei 

Likewise one also derives the width due to the i — 2 term as 

f ^ 1.9 X 10- . (.*/3 X X { |l _ ^3^3,,,_ I > l (32) 

The total width 59/9 at a given 9 is then the quadrature sum of the contributions (at that 
value of 9) from both the i — 1 and ^ = 2 terms, as in Eq. (21). 

A graph of 5^/^ = 59/9 versus I = n /9 for the i = 1, i = 2 components and their total 
is shown in Figure 5. We explained in the end of section 3 that the smoothing kernel for 
59/9 is Si one-sided gaussian function in the case of time delay perturbations, i.e. structures 
on the LSS are invariably biased towards having smaller perceived angular sizes irrespective 
of how time delay may tilt the LSS. Specifically the kernel K{9, £) that enters the integrand 
for the lensed CMB temperature correlation function, viz. Eq. (A5) of Seljak (1996) is 

M».Q = {„_ ',4; (33) 

where 0.6a{9) = 59 with 59 given by Eqs (31) and (32) added in quadrature, and the factor 
of 0.6 is due to the fact that the standard deviation of a one-sided gaussian having the usual 
form for its exponent has standard deviation 0.6a. 

Thus, after applying the (normal) simplifying procedures of Limber approximation and 
retention of only the isotropic term Jo{£9) in Seljak's (A5), one arrives at the lensed corre- 
lation function due to time delay 

Ce= 9d9 / idiK{9,i)CiJo{i9), (34) 
Jo Jo 

which is to be contrasted with Eq. (A7) of Seljak (1996). Owing to the smallness of the 
effect as dictated by the numerical coefficients in Eqs. (31) and (32), however, we did not 
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convolve the acoustic peaks according to the recipe of Eq. (34), because the distortion is 
completely unnoticeabl^. We can therefore corroborate HCOl, after this comparatively more 
detailed and specific investigation (see section 1), by concluding that relative time delay by 
foreground primordial density fluctuations does not lead to any appreciable changes in the 
prediction of the standard cosmological model on the CMB acoustic anisotropy. 



^In order for any smearing effects to distort the acoustic peaks to the point that the resulting standard 
cosmological model no longer matches WMAP3 data, it is necessary for <j{9)/9 to be at least 10 % at the 
center of the first peak. 



100 200 300 400 500 600 700 

Multipole moment i 

Fig. 5. — The percentage variation of the size of an acoustic harmonic i due to time delay 
differences along the LSS, i.e. Eqs. (31) and (32), is plotted as a function of £ for contribution 
from the mode i = 1 (long dashes) and i = 2 (short dashes) of primordial density fluctuation 
in Eq. (1) and Table 1. The total 5^/^ curve, obtained by adding the two dashed lines in 
quadrature, is shown as the solid line. Note that the i = 3 mode is ignored because of the 
smallness of its effect. 
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7. Can matter clumping ELfFect the 'mean' geometry of space? 

In spite of the feeble acoustic distortion induced by the time delay perturbation of the 
primordial matter distribution, its effect of biasing the acoustic peaks towards higher values 
of £ (as compared with the values under the scenario of a homogeneous Universe) remains 
in principle an interesting phenomenon. To illustrate, let us repeat the acoustic smearing 
computation in the previous section using a kernel K{9,i) as determined by assuming a 
larger normahzation of the matter power spectrum (than the WMAP value of Eq. (3)), viz. 

5$ fti 0.057. (35) 

The consequence, shown in Figure 6, is no longer neghgible as before. The systematic shift 
of the acoustic peaks towards higher i (or smaller 9) is now evident, and tends to force the 
best-fit model parameters to assume new values - that involve Q < 1 in particular. Thus 
clumping does not just introduce local fluctuations in the metric tensor (as first noted in 
section 3 of Einstein 1917), they can mimic space curvature, by altering the mean geometry 
of space as revealed by the statistical behavior of remotely emitted light along many and 
varied directions. The magnitude of this influence is by no means necessarily small: even 
the value of 5$ in Eq. (35) lies well within the perturbation criterion of 5$ <^ 1, yet it led 
to a finite CMB 69/ 9 which, for the i — 1 term of Eq. (1), is not <^ 1 because it maxmizes 
at 32.5 % when Eq. (31) is rescaled on the right side according to 69/9 ~ (5$)^. 

In the light of the above development, let us return to the question of why global space 
as measured by WMAP has zero curvature. It has widely been accepted that inflation 
solved the flatness problem, but does not answer why 6^ is as small as the value in Eq. (3). 
which implies an unusually weak scalar-field coupling. This view is still correct, so long as 
the 'flatness' is understood as an intrinsic rather than observed property of space, because 
indeed 5$, no matter what values it assumes, does not alter the mean intrinsic curvature. 
But since, as demonstrated in this work, even within the perturbation regime a larger (5$ 
can make an intrinsically flat geometry appear to WMAP as, on average, curved, then one 
must acknowledge that inflation has not explained why the observed mean curvature is zero. 
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Fig. 6. — Under the assumption of an enhanced spectral normalization 5$ for the primordial 
density perturbation, i.e. with Eq. (35) replacing Eq. (3), the rescaled (hence mock) WMAP3 
data and rescaled WMAP3 standard model prediction of the acoustic peaks are plotted, with 
the latter as a dotted line. This model is then re-processed by the lensed correlation function 
due to time delay, viz. Eq. (34) of section 5, again with the new 5$ value of Eq. (35). 
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A. APPENDIX 

The basic relationship between gravitational potential perturbation at veirious 
length scales and the matter power spectrum 

During some epoch t, or redshift z, let the typical random excursion of the gravitational 
potential (from its mean value in a given cosmology) be at wavenumber k or lengthscale 
2Tra/k. By the Poisson equation, is expressible in terms of the matter density variation 
Sk{z) = Sp^{z)/pm{z) as = A'KGpm{z)a^5k{z)/k'^, where a = a{t) = 1/(1 + z) is the 
expansion factor. Squaring both sides, we obtain the variance the additive nature of 
which enables us to write the equation in differential form, as 

d^l = [ATTGpUzWa'^. (Al) 

Under a matter dominated scenario, Sk oc {1 + z)~^ ex a, and since Pm{z) = (1 + z)^pm{0) — 
Pm{^)/ci^-i we see that is independent of epoch. Thus it may be written in terms of the 
present day quantities 

d^l = [47rGp^(0)]^^, (A2) 

where d5l{Q) relates to the z — Q matter power spectrum P{k) as inferred from WMAP's 
standard model parameters (or, better still, the WMAPl/2dFGRS external correlation, see 
below) by 

,Sli0) = f'^. (A3) 

It should be noted, however, that if vaccum energy (or dark energy) rather than matter 
dominates the low z Universe, Eq. (A-2) would have slightly overestimated the value of $fc 
at an earlier epoch. 

Before proceeding further, therefore, we must estimate a strict limit on the inaccuracy 
of Eq. (A-2). Under a A-dominated scenario, the rapid expansion freezes the density distri- 
bution - the masses cannot move fast enough to counter it. As a result, 6k tends to a constant 
independent of a (i.e. no further growth of density contrasts in the linear regime), and $fc 
decreases as a~^. In a Qjn — 0.3, Qa — 0.7, and h — 0.7 cosmology (Bennett et al 2003, 
Spergel et al 2007), vacuum domination occurs at about z 0.32 (when ^1^(1 + z)^ — JIa). 
Thus, between z ^ 0.32 and z — 0, d^k decreased from its hitherto constant value by the 
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fraction 1 — a ^ 30 % (where a refers to the expansion factor at 2; = 0.32). It is clear then 
that any hghtpath integrations of the potential from z = back to some remote past (or 
vice versa) will certainly not lead to an overestimation of by more than 20 % if ones uses 
Eq. (A-2) for <^k- 

Let us for the time being suppose that the z — matter power spectrum has the simple 

form 

P{k) = Ake-^'', (A4) 

so that at small k the power spectrum takes the Harrison-Zel'dovich form P{k) ~ k. Then, 
assembling Eqs. (A-2) and (A-3), we have 

d^l = SG^plAe-'"^ = {6^)^Ae-"'^, (A5) 

where 

(5$)2 ^ M = ^^^^A, (A6) 

with Hq being the Hubble constant and c = 1 here-and-after. At sufficiently large k, cor- 
responding to wavelengths smaller than the size of the Universe (in today's distance scale) 
during matter-radiation equipartition, P{k) cuts off because the modes could not grow. 
Thus, in the simple manner by which Eq. (A-4) depicted the matter spectrum, 6 ~ 33 Mpc, 
the equipartition horizon for a fiat Universe with h = 0.7. The advantage of Eq. (A-4), 
as can be seen in sections 3 and 4, is that it affords us an analytical treatment of CMB 
foreground re-processing by primordial matter - time delay and lensing in particular - by a 
method of successive approximation from low to high orders which reveals unamiguously the 
coherent spatial scales of these effects. 



B. APPENDIX 

Angular correlation function of deflection due to primordial foreground matter 

The correlation function C{6) as defined by the first equality of Eq. (25), section 4, 
may be computed using Eqs. (24) and (11). We find 



X 



- 2x^9^ - + {3x^9^ - 2r2) J arctan (f) 

^ (Bl) 

where x — {x' + x")/2, x — x' — x" as in section 3, and 

r^^x^ + x^^^. (B2) 
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because with the new x- variables we have, from Eq. (12), r — (x, x6). By setting ^ = one 
arrives at the variance of the deflection angle for a single light ray. 



Co = C(0) = 



,2 ex 



dx 



dx 



[x — xy — jxP' 



x" 



2x^ — 2hx arctan ( ^ 



/ ^2 



1 arctan | — 

X 



(B3) 



Evaluating the integral yields 



(B4) 



in the limit x^h. 



In general, the cross-correlation function C{d) of Eq. (B-1) may be expanded as a Taylor 



series: 



(B5) 



where 



Co = 



32(5$ 



2 rx 



dx 



JO 



' dx^^^—^ 4 

^ x^ 



- 1^21^2 r 



b ^ 62 
6 — arctan 2 — 



X 



X 



Z.2 



X 



+ 62 



I^io 



3— arctan — — 2 — -r^^ — -r^ 

X X x^ + 



(B6) 



Again as before, in the limit x ^ h only the leading term needs to be kept. The result is 



(B7) 



We finally arrive at the expansion 



x2 3 6 



1 X2g2 

20 62 ^ 



(B8) 



after substituting Eqs. (B-4) and (B-6) into Eq. (B-5). 
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